HIGHER-ORDER LAPLACE EQUATIONS 
AND HYPER-CAUCHY DISTRIBUTIONS 



ENZO ORSINGHER AND MIRKO D'OVIDIO 



Abstract. In this paper we introduce new distributions which are solutions 
of higher-order Laplace equations. It is proved that their densities can be 
obtained by folding and symmetrizing Cauchy distributions. Another class of 
probability laws related to higher-order Laplace equations is obtained by com- 
posing pseudo-processes with positively-skewed Cauchy distributions which 
produce asymmetric Cauchy densities in the odd-order case. A special at- 
tention is devoted to the third-order Laplace equation where the connection 
between the Cauchy distribution and the Airy functions is obtained and ana- 
lyzed. 



1. Introduction 

The Cauchy density 



pM= \^tw) (L1) 



solves the Laplace equation (see Nane [5]) 



dt 2+g^ = °> xeR,t>0. (1.2) 



The n-dimensional counterpart ol ( |1.1[ ) 



p(x,t) = — ^ s-, xeM^.OO (1.3) 

7T* ( t 2 + | x |2)2 

with characteristic function 

/ e i ( a '*>p(x,t)dx = exp(-t\a\) (1.4) 
it"- 1 

solves the n-dimensional Laplace equation 

j'=i 3 

The inspiring idea of this paper is to investigate the class of distributions which 
satisfy the higher-order Laplace equations of the form 

^7 + ^-0, ,6l,t>0 (1.6) 
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In a previous paper of ours we have shown that the law 



Pi(x,t) 

solves the fourth-order Laplace equation 



t x 2 + t 2 

+ 1 4 



(1.7) 



d 4 u d A u 



t>0. (1.8) 

In Section 2 we analyze distributions related to equations of the form 

9 2 "m d 2 "u 







(1.9) 



dt 2 " dx 2 " 

which can be expressed in many alternative forms. The decoupling of the 2™-th 
order differential operator in (1.9) 



dt 2 " 



d 2 " 
dx 2n 



2™" 



n 



fc=_(2 n - 1 -l) 
k odd 



dt 2 



+ e 2" 



d 2 
dx 2 



suggests to represent distributions related to (1.9) as 



1 



E 



- 2 "-\ = _ ( ^- 1 _ 1) * 2 +(^r 

k odd 



n > 2. 



(1.10) 



that is the superposition of Cauchy densities at imaginary times. Alternatively, we 
rive a real 



give a real- valued expression for (1.10) as 
t{x 2 +t 2 ) 



2 »-i_x 2"~-l 



2 Tr(x 2n +t 2n ) 



k—1 j odd 

k odd j^k 



E II (x 4 + t 4 + 2x 2 t 2 cos 



(1.11) 



The density ( |1.11[ ) can also be represented as 
t(x 2 +t 2 ) r ^- 1 



71 "l X 

k odd 



'ifi 2x 2 t 2 cos 



n > 2. 



(1.12) 



Each component of the distribution ( 1.12 ) is produced by folding and symmetrizing 
the density of the r.v. 



V(t) = C tcos 



kn 
2™ 



t sin 



2" ' 



<> 0. 1 < A; < 2 



n-l 



1, odd 



where C(i), t > is the Cauchy symmetric process. The distributions (1.12) differ 
from the Cauchy laws since they have a bimodal structure for all n > 2 as figures 
below show. For n = 2, the distribution (1.11) reduces to ( |1.7[ ) if we assume that 
the inner product appearing in formula ( 1.11 ) is equal to one. Of course, the density 
( |1.12[ ) coincides with ( [L~7| ) for n = 2. For rt = 3 we get from ( |1.11[ ) and ( |1.12[ ) that 

t(x 2 +i 2 ) 



P2z(x,t) 



'V^Tr(x s +t») 

ttx 2 + t 2 ) r 



(x 4 + t 4 - V2x 2 t 2 ) cos J + (x 4 + t 4 + V2x 2 t 2 ) sin J 
8 8 



cos 
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In Orsingher and D'Ovidio [TT] we have shown that the density (1.7 1 is the proba- 
bility distribution of 

Q(t) = F{T t ), t > 

where F is the Fresnel pseudo-process described in [TT] and T t , t > is the first 
passage time of a Brownian motion independent from F. We note that 

Q(t) = F(\B(t)\), t>0 
has density coinciding with the fundamental solution of the fourth-order heat equa- 
tion 

du d A u 
di = ~dx±~ 

We prove also that, for k G N, there are non-centered Cauchy distributions which 
solve the equations 



0. 



Qfik+1 dx 2k + l 

If X2k+i{t), t > is the pseudo-process whose density measure 

H2k+\{dx,t) = n{X 2 k+i{t) G dx} 
solves the heat-type equations 



(1.14) 



du 

at 



Q2k+1 



k e N 



Q x 2k+1 ' 

and S_l_ (t), t > is a positively skewed stable process of order 2/ ; ( 

Pr{X 2k+1 (S^_(t)) g dx}/dx = — tcos W^T) 

(x+ (-l) fc+1 tsin 2(2 fc +1) j +t 2 cos 2 



(1.15) 
we have that 



2(2fe+l) 

(1.16) 

We show below that the densities ( 1.16 1 solve also the following second-order p.d.e. 

d 2 u 



d 2 u d 2 u 



2 sin ■ 



dt 2 dx 2 2(2k + l)dtdx' 

We have investigated in detail the case of third-order Laplace equation 

d 3 u d 3 u 
dt 3 dx 3 

and have shown that 







Pr{X 3 (Si(t)) G dx} =dx / - TJ =Ai 
3 Jo v3s 



''is 



t 1 



, y/3 X-t y/3 

=dx 1— = dx 



/3s 
t 



-M 



t 



/3s 



ds 



(1.17) 
(1.18) 



=dx 



-t 3 
t cos -| 

6 



2 x 2 + xt + t 2 



(.x + isinf) 2 +t 2 cos 2 f 



The pictures of the Cauchy distributions (1.16) show that the location parame- 
ter t sin ; 



t cos 



2(2fc+l) 



2(2fc+l) 



tends to zero as k — > oo while the scale parameter tends to one, 
t. This means that the asymmetry of the Cauchy densities decreases 
as k increases and is maximal for k — 1. The decrease of parameters of (1.161 (with 
k increasing) is due to the growing symmetrization of the fundamental solutions of 
equations (1.151. 
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By suitably combining the distribution (1.16) for k = 1, we arrive at the density 

(x 2 + t 2 ) cos I + xt 



P6{X,t) = —r-t 



2 2 tt (x 2 +t 2 + irfcos§) 2 + 2x 2 i 2 cos§ 



which solves the equation 



d 6 u d 6 u 
dW + dx^ 



= 0. 



(1.19) 



(1.20) 



The probability density (1.19) displays the unimodal structure of the Cauchy dis- 
tribution. 



2. Hyper Cauchy distributions 

In this section we analyze the distribution related to Laplace-type equations of 
the form 

a 2 " <9 2 " \ 

For n > 2 we obtain a new class of distributions having the form 

t{x 2 +t 2 ) 



(2.2) 



where g(x,t) is a polynomial of order 2™ — 2 2 . For n — 2, formula (2.2) yields the 
distribution 

t(x 2 + t 2 ) 



Pi(x,t) = 



x e R, t > 



(2.3) 



V2ir(x 4 + t 4 ) ' 

emerging in the analysis of Fresnel pseudo-processes (see Orsingher and D'Ovidio 

ELI]). 

The main result of this section is given in the next theorem. 
Theorem 2.1. The hyper Cauchy density 



P'2 I'-') = 7r 2^— T E 



k=-(2 

k odd 



te 1 ^ 
x x 2 + (te 1 ^) 2 



solves the equation 



d 2 



if 



dt 2 " dx 2 " 



u = 0, x e K, t > 0, n > 1. 



(2.4) 



(2.5) 



A real-valued expression of (|2.4|) reads 
t{x 2 + t 2 ) 



E 

/c=i 

odd 



A'TT 



2 n- 



COS ■ 



k^j = l 
j odd 



[] (x 4 +i 4 + 2^ 2 cos-^ T 



or equivalently 



P2"{x,t) 



t{x 2 +t 2 ) 



E 



cos 



2™- 2 7T ^ x 4 +t 4 + 2x 2 t 2 COS ^r-r 

k odd 



(2.6) 



, for n > 1. (2.7) 
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Proof. In order to check that (2.4) satisfies equation (2.5) wc resort to Fourier 
transforms 

f + oo 



e l0x u{x,t)dx. 



Equation (2.5) becomes 

d 2 "U 



dt 2 



+ {-i(if n U 



d 2 "U 
dt 2n 



+ /3 2 U = 0. 



The solutions of the algebraic equation associated to (2.8) have the form 
r s = |/3|e l7r ^, 0<j<2 n -l. 



In order to construct bounded solutions to (2.8) we restrict ourselves to 

2" 



1709, t) 



1 



L -i 



in-l X! e 
fc = -(2"- 1 -l) 
k odd 



-t\P\e2™ 



(2.8) 
(2.9) 

(2.10) 



where the normalizing constant in ( 2.10 ) is chosen equal to 1 /2 n ~ 1 so that U (/3, 0) = 
1. The inverse of (2.10) is (2.4). We check directly that each term of (2.4) has 
Fourier transform solving equation (2.8). For all odd values of k, we have that 

10X 



„i/3z I ° 



if 



\dt 2 " da 



x 2 + (te 1 ^) 2 



dx 



& 



2" 



2™ o2 T 



dt 2 

= (/3 2 V fe7r -M 
= ((-l) k f3 2 " +/3 2 



-t\P\e 



-t\f3\e'2^ = Q 



because k is odd. In order to obtain (2.6) we observe that, in view of (2.4) we can 
write 



1 



p 2 « (&,*)=- Y] 

7T ■ ■* 



Ck 



^|2fe-l| a; 2"-|2fe-l|-l 



k odd 



k odd 



where 



TJ (x 2 + (te l ^ 2 ) = x 2 ' 1 +t 2 



(2.11) 



/fe=-(2 n-1 -l) 
k odd 



and Cfe are constants evaluated below. Result (2.11) can be obtained directly by 
solving the equation x 2 + t 2 = or by successively regrouping the terms of the 



right-hand side of (2.11). We have at first that 

2 



L -i 



L -i 



n (x 2 +(te^f)= n 



k— — (2™ — — 1) 
k odd 



k — 1. k odd 



t 4 + 2x 2 t 2 cos 
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fe=l, k odd 



x 2 " + t 2 " + 2x z i z cos 



2 + 2 



=x 2 " + e 



In view of (2.11 ) we can rewrite (2.4) as 

2 n-l_ 1 2" _1 -1 



p 2 r>(x,t) 



t 



n2 n - 1 (x 2 " +t 2n ) 



fe= _(2«-i_l) :! -=-(2»- 1 -l) 
fc odd j odd,j^/c 



where 

2 »-i_ 1 2" _1 -1 



fc=-(2"- 1 -l)j = -(2"- 1 -l) 
fc odd j odd, j^k 

2 n ~ 1 -l 2 T1_1 -1 



J2 II (^ 4 +^ + 2^ 2 cos 2 ^ T )(^ + (te-^))e l 



fc=-(2"- 1 -l) 3=1 

fe odd j odd, jjtk 

2 «-i_ 1 2" _1 -1 



E H (a; 4 + i 4 + 2x 2 t 2 cos^ T )(^e^+i 2 e-^+x 2 e-^ + ^e^) 



fc=l 3=1 
fc odd j odd, jj=k 



kn 



-2 { x 2 +t 2 ) y: cos^ n 



(.T 4 + t 4 + 2l 2 t a C0B-^ T ) 
v On— 1 y 



fe=l 
fc odd 



3=1 
j odd, j^k 



and thus 



t(x 2 + t 2 ) 



kit 



• - y cos^ TT (:r 4 + i 4 + 2x 2 i 2 cos-^ T ). 



A: odd 

Furthermore, from the fact that 

2 n-2_ 1 



3=1 
3 odd, j^fc 



X 2 " + i 2 " 



fc— 1, k odd 



n U 8 +i 8 + 2^ 4 cos^ 



we obtain that 



t(x 2 + 1 2 ) 

2™- 2 7T 



E 

fe=i 

fe odd 



kiv 



c 4 + t 4 + 2x 2 t 2 cos 



□ 
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Remark 2.2. In order to prove that the density ( 2.6 ) integrates to unity we present 
the following calculation 

,+oo x 2 +t 2 ,+oo 

/ ; — ax —2 

J-oo x i + t 4 + 2x 2 t 2 ' 1 



x 2 +t 2 

; — ax 

2x 2 t 2 cos 

y 2 + i d 

!/ 4 + l + 2ycoB 5 £ r V 
1 



COS 



X 4 +t 4 



2 



f DO 



d0 



tan 4 0+1 + 2 tan 2 cos ^ cos 4 



2 ft d0 
* Jo sin 4 6» + cos 4 9 + 2 sin 



1— cos — 

1 STL 



2 sin 2 (9 cos 2 9 cos ^ 

rZ0 



2^ sin 2 20 



d0 



1 f 2 



l Q 1 _i (l _ cos ^ T)( i^) 



./() 



d(j) 

1-COS -^r , , 

1 i^T 1 



4^ + \ (l -COS COS c 



2t 

2 ^ (/(.•■■ 

~* io (3 + cos ^fer) + (1 - cos 2^-) COS ( 
2 2tt 

* ^(3 + cos ^t) 2 - (1 - COS 
_ny/2 1 



t 

IT 1 



/!• 



- COS 



in 



* cos ' 

view of (2.12), we can conclude that 

p + OO 



(2.12) 



From ( p?7| , 

Remark 2.3. From ( |2.4| , for 

, \ 1 

MM) = — 

with Fourier transform 



/+oo 
p 2 ™(x,t) dx = 1 
-oo 

n = 2 we obtain that 

1 



te 1 * fe l <i 

.t 2 + (fe*f) 2 x 2 + (<e~ i f) 



From (2.6) and 
The law 



( pj| we 



f + e^ x p 4 (M)<fe = e~^ l/3| cos4^ 
-oo V 2 

have that 



, . t x 2 + t 2 
(2.131 has two maxima as Figure [l] shows. 



(2.13) 
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Remark 2.4. For n = 3, from (2.6 1, we have that 
t(x 2 +t 2 ) 

From the fact that 



(x 4 + t 4 + 2x 2 t 2 cos cos 



3tt 



i 4 + 2x 2 t 2 cos ^ ) cos 



cos 



37T 7T 37T 7T 

= — cos — and cos — = sm — 
4 8 8 



we write 



i(x,t) 



t(x 2 +t 2 ) 
Z 2ir(x 8 +t & ) 



(x 4 + t 4 + V2x 2 t 2 ^) sin | + (a; 4 + i 4 - \/2izi 2 t 2 ) cos | 



From (2.7) we have also that 

x 2 +t 2 



p s (x,t) 



t 

'2k 



x + t 2 



c 4 + t 4 - V2x 2 t 2 



t 4 + y/2x 2 t 2 8 



(2.14) 



From (|2.4|) we obtain the characteristic function 

(^ S in|) 



e tl3x p 8 (x,t)dx 



2 2 



e -t|/3|c OS | C0SI 



-t\p\ sin ■ 



cos ( tp cos 



The density p% (x , i ) is a bimodal curve as well as p^ (x , t) . The maxima of p% (x , t) are 
heigher than those of pa(x, t) as Figure [I] shows. Also P2™{%, t) displays a bimodal 
structure with the height peaks increasing as n increases. The form of P2 71 i x , t) 
reminds the structure of densities of fractional diffusions governed by equations 



d v u X 2^ 2u 
dx 2 



for 1< v < 2 (see [ID]). 



Remark 2.5. The result (|2.7|) can conveniently be rewritten as 

x 4 + t 4 + 2xH 2 



P2" (X,t) 



t 



ir(x 2 +t 2 ) 



— y 

n-l 



kn 



k=l 
k odd 



t 4 + 2x 2 t 2 « 



—, — cos — 

kit 2 r 



(2.15) 



The factor in square parenthesis measures, in some sense, the disturbance of p2" 
on the classical Cauchy. For n — 2, we have in particular that 



p 2 2 (x,t) 



t 



1 



7T(X 2 +t 2 ) y/2 



1 



2x 2 t 2 



x 4 + t 4 



t x 2 +t 2 

7^x 4 + t 4 ' 



(2.16) 



The density (2.161 has two symmetric maxima at x = ±tv \/2 ~ 1 and a minimum 
3rsinghe 

9k(x,t) 



at x = (see Fig. 6 of Orsingher and D'Ovidio [H]). The terms 

x 4 + t 4 + 2x 2 t 2 



x 4 + t 4 + 2x 2 t 2 cos 



(2.17) 



display two maxima at x = ±t with height depending on k and whose profile is 
depicted in Figure [2] 



HYPER-CAUCHY LAWS 



Figure 1. The profile of the functions P4 (dotted line), formula 
(|2.13[) and p 8 , formula (|2.14[). 
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Figure 3. The profile of the functions P2», formula (|2.7[), for n = 5, 10, 15,20. 



Remark 2.6. The density p2^(x,t) can be written as 

t{x 2 +t 2 ) 



P2" (X,t) 



Q(x,t) 



(2.18) 



2 n - 2 Tr(x 2 " + t 2 ") 

where Q(x,t) is a polynomial of order 2™ — 2 2 . For n — 2 the function Q(x,t) 
reduces to cos j. For n = 3, 

Q(z, t) =(.x 4 + t 4 + V2x 2 t 2 ) sin J + (x 4 + t 4 - V2x 2 t 2 ) cos J. 

8 8 



The expression (2.18) shows that the probability law p 2 n(x,t), x <E R, t > shares 
with the classical Cauchy density the property of non-existence of the mean value. 



Remark 2.7. The density of the hyper Cauchy can also be presented in an alter- 



native form by regrouping the terms in the right-hand side of (2.7) as 

2 n-2_ 



E 

k=l 
k odd 

E 

fc=i 

k odd 



sin 



COS 



x 4 +t 4 - 2xH 2 cos gfer x 4 + t 4 + 2x 2 t 2 cos ^ 



(2.19) 



1 (x 4 + t 4 + 2x 2 t 2 cos 2^-) sin ^ + (x 4 + t 4 - 2xH 2 cos ^^r) cos 



c 8 + 1 8 - 2x 4 i 4 cos 
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For n — 3, from (2.19), we get again that 



t(x 2 + t 2 ) 



2n(x 8 + t s ) 
Remark 2.8. The r.v 



(V + t 4 + V2x 2 t 2 ^) sin ^ + (V + t 4 - v^ 2 ) cos 



fc7T 



W(t) 



7rfc\ . 7rfc 
G t cos — — i sin — — 
V 2" / 2™ 



(where C(t), t > is the Cauchy process) has probability density 

. 2t{w 2 + t 2 ) cos 

fk{w,t) = —— — — — — — j——, w>0. 

ir(w 4 + t 4 + 2w 2 t 2 cos ) 



Indeed, we have that 

Pr{W{t) < w} 



-\-w-\-t sin ^ 



dy 



t cos 



—iu-\-t sin £ 



% 7r(?/ 2 +t 2 cos 2 |f) 



(2.20) 



(2.21) 



(2.22) 



and 

fk(w,t) 



dw 



-Pr 



C t cos 



Trfe 



f sin ■ 



71 A - 
2" 



< w 



t cos 



fe7T 



tcos^ 



7T Uw + t sin ^) 2 + t 2 cos 2 7T ^(-io + * sin §£) 2 + t 2 cos 2 §f J 



i cos 



t COS 



fc7T 



7T u; 2 + 2u;tsinP+t 2 tt [ w 2 - 2wt sin + t 2 



2t(w 2 + t 2 ) cos |S 



~7r(u; 2 + t 2 + 2wtsin|^)(w; 2 + i 2 - 2wtsin^) 



2i(ui 2 + i 2 ) cos 



7r(u) 4 + t 4 + 2w 2 i 2 cos^fer) 



because 



2 sin 2 



kit 
2™ 



1 — cos ■ 



k tt 

in— 1 ' 



By symmetrizing (2.20) as follows 

Z(t) 



Wi(t) - W 2 (t) 



where Wi(t), W2(i) are independent copies of W(t) we obtain a distribution of the 
form 



hk(w,t) 



t(w 2 + t 2 ) cos |£ 



ti(u; 4 + i 4 + 2u; 2 i 2 cos ^f^; 



(2.23) 



which coincides with each term of (2.15). This construction explains the reason for 



which each term in (2.15) has two symmetric maxima at w = ±iy 2 sin |£ — 1 for 
k : sin f| > \. 
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Figure 4. The figure shows how the distribution (2.231 can be 



constructed from the Cauchy density by folding and symmetrizing, 
in the cases n = 3, k = 1 (top figures) and k = 3 (bottom figures). 



The dotted line gives the density of the folded distribution (3.13) 







3. Higher-order Laplace- type equation 
Let us consider the pseudo-processes related to higher-order heat-type equations 



x G 



t > 0, n > 2. 



(3.1) 



du d n u 

where c„ = (— 1)5 +1 for n even and c„ = ±1 for n odd. 

Pseudo-processes constructed by exploiting the sign-varying measures obtained 



as fundamental solutions to (3.1) have been examined in many papers since the 



beginning of the Sixties. A description of the procedure of construction of pseudo- 
processes can be found, for example in Krylov 4J, Ladokhin 6J, Hochberg [3], 
Orsingher [S], Lachal [5]. In the case where n = 2k+l, c 2 fe+i = — 1, the fundamental 



solution to (3.1) reads 

u 2k +i(x,t) 



1 

27 



+00 



i -i l3x +i(-l) k tl3 2h + 1 d o 



In particular, for k = 1 



u 3 (x,t) 



1 



cos (#e + f3 3 t) d@ 



-.Ai 



x 

Wt 



(3.2) 



(3.3) 



where 



Ai(x) 



-x 
3 



3/2 
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is the Airy function (see for example Lebedev [7]). 

In this section we study the composition of pseudo-processes with stable pro- 
cesses S a (t), t > 0, a € (0, 1) whose characteristic function reads 



Ee */3S Q (t) = exp ^_ t ^a e -i^^ = exp f- at \f3\ a ( 



exp ( -at\p\ a ( 1 -i9^- tan ^ 



(3.4) 



where a — cos 717/2 > and 



cot(-)tan(^) 



The parameter 7 must be chosen in such a way that 9 £ [—1, 1] for a € (0, 1). The 
skewness parameter # = 1 (that is 7 = a) corresponds to positively skewed stable 
distributions. For the density 



1 

2^ 



e -ipx Ee ipS a (t) d p 



(3.5) 



we have the scaling property 

For a £ (0, 1), we have the series representation of stable density (see [101 page 
245]) 

Pa (x; 7 , 1) = I EC-ir ^^ 1 ^ ^"^' 1 sin (f (7 + «)(»•+ 1)) ■ ( 3 - 6 ) 

r=0 

Theorem 3.1. TTie composition of the pseudo-process X2k+i(t), t > wit/i i/ie 
stable process S i _ (t), t > 0, k £ N, /ias a Cauchy probability distribution which 
can be written as 



Pr{X 2k+1 (S^_(t) £ dx}/dx 



t cos 



2(2fc+l) 



2(2fc+l 

(3-7) 

with x £ K, t > 0. TTie density function (3.7) is a solution to the higher-order 
Laplace equation 

W^T + dx^I =0 > x ^ R > t>0 ( 3 - 8 ) 
Proof. For = 1, a = 7 = l/2k + 1, in view of (3.4) we have that 

/+00 
e^Pr{X 2k+l {S^_{t))£dx} 
-OO 

poo />+oo 

= / Pr{S 1 (t) £ ds} / e i/3:l; U2fc + i(x,s)da; 
Jo 2k+1 J-00 



Pr{Si(t)£ds} 



exp —t 



k o2fc+l 



2fe+l 7T 
COS 



exp — i|/3| cos 



2(2fc + l) 



2(2fc + 1) 
i(-l) fe T ^ T sin 



1 — i si 



m ((-l)^ 2 ^ 1 ) 



tan 



2(2fc + 1) 



2(2fc + 1) 
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, ^'*-^ |COS 2(2^- j( - 1) ^ Sm 2(2fcTr) 



(3.9) 



This is the characteristic function of a Cauchy distribution with scale parameter 
t cos 2 (2fc+i) an d location parameter t(— l) k+1 sin 2 (2k+i) • Formula (3.9) can also 
be rewritten as 



U(p,t)=exp \-t\, 
= exp ( -t[ 



cos 



COS 



exp (-t\p\e l 2 ( 2fc+i)( i)* 1 



2(2fc + l) 

7T 

2(2£; + 1 
§_ 



sm 



\P\ 2(2* + 1) 



— z sm 



2(2fc + l) 



(-1) 



k /3 



The Fourier transform of equation (3.8) becomes 

FP k+1 T7 

i j i + (_, /? )2fc+l, 7 = 0. 



Q t 2k+1 



(3.10) 



(3.11) 



The derivative of order 2k + 1 of (3.10) is 



f)2k+l n 



g 2(2fc + 



2k+l 



U(j3,t) (3.12) 



and this shows that the Cauchy distribution (3.7) solves the higher-order Laplace 
equation (3.8 1 . □ 



Remark 3.2. We notice that 

rOO i poo 

/ Pr{X 2k+1 (S i (t))edar}=- / 

2 V 2fc + 1 



1 + 2/ 2 



(3.13) 



which is somehow in accord with Lachal [5]. The results (3.7) and (3.13) show that 
the mode of the Cauchy law (3.7) approaches the origin as k increases. 



Let us consider the process of the form X 3 (Si (t)), t > where X 3 is a pseudo- 
process whose measure density is governed by the third-order heat equation 

s— " R ' i>0 < 3 - 14 > 

and Si is the stable process of order 1/3. The distribution of X 3 (Si (£)), t > 
reads 



f°° 1 

Pr{X 3 (Si(t)) £ dx} = dx / ^=Ai 
3 Jo v3s 



t 1 



where 



for which 



Pr{Si ft) G ds} = ds-^=Ai 
3 s v/37 



/3s 



v 3s 

s > 0, £ > 



ds (3.15) 
(3.16) 



/•OO 

=(w = i/v 7 ^) = 3 / Ai(iu) du; = 1. 
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Corollary 3.3. The law (3.151 solves the higher-order Laplace equation 

d 3 u d 3 u 



dt 3 dx 3 



o, x g R, t > o 



and can 6e written as 



dx 



Pr{X 3 (Si_(t)) G dx} = 

3 7T 



\/3 . 

2 11 

n 2 



3*1 
4 



dx 3 1 / 2 



7T 2 X 2 + xt + t 2 

3 1 / 2 ; x-t 



-dx 



2tt x 3 - t 3 ' 



(3.17) 



(3.18) 



Proof. The Fourier transform of (3.15) becomes 



e i^p r {x 3 (Sx(t)) edx}= I , 

oo 3 JO 



Pr{Si(t) G ds}. 



(3.19) 



We show that (3.16) is a stable law of order 1/3. In view of the representation of 

(3.20) 



the the Airy function (4.10) of Orsingher and Beghin [10 

27 

~3 



3-2/3 ~ (3l/3 w )* . / 27r \ /jfe + 1 

A *H = — — E fcj sm ( — + ) r 



fc=0 



we can write that 



t 1 



s V3s 



/3s 



t 

3irs\/s 



E 

fc=0 



/ \ fc 1 
fcl 



27T„ A /fc+1 

Mil i — (fc + 1) r 



We consider the series expansion (3.6) of the stable density (with t — 1) for which 
(JO) holds true. For a = 7 = 1/3 (that is = +1), x = s '/t 3 in EM we get that 



s 1 

"i V^ ; 3' 



fc=0 

=(by 4.5 of [TO] 
1 t 4 ' 



' -•') 3.E (.,■,) 



sm 



fc + 1 



E 



37T S-c/s » - 3/ 

V fc=0 



1 /27T,, ,\ /fc+1 

sin [ —(fc + i) r 



fc! 



=f 1 



t 1 



/3s 



and thus, from (3.5|, we have that 
1 



t 3Pl * \t 3 ' 3 



s 1 \ / 1 \ i 1 . 



s V3s 



'3s 



s,t > 0. 



We now evaluate the integral (3.19). We have that 



10 
= exp 



e~* p s Pr{St(t) G ds} 

(-cob Jt| -jS 3 ]* (l-isgn (-/? 3 )tan|)) 
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= exp 



V3. 



t\P\ 1 + zsgn 08) 



73 



: exp L^t^l-i^ 



(3.21) 



since sgn (— j3 3 ) = sgn (— /?) = — sgn (/3) = — JL. From (3.211 we infer that 



dx r + °° 

Pr{X 3 (S h (t))ed X }=-j^ 



e-^expf-^^l-^j d/3 



dx 



dx 3 1 / 2 



t 



if if 7T 2 a; 2 +a;i + t 2 



[X-, 2 j i 4 

3 1 / 2 * x-t 
2tt - t 3 



(3.22) 



□ 



Remark 3.4. We observe that the r.v. X^(Si(t)) possesses Cauchy distribution 

with scale parameter \/3t/2 and location parameter —t/2. Furthermore, it solves 
the third-order Laplace-type equation 

d 3 u d 3 u 



dt 3 dx 3 



0. 



Remark 3.5. From the fact that 



1 



x 

3^3* 



1 fx 



( 2 x 3 ' 2 



3vr V t 1/3 \3 3 / 2 y/t 



, x,t > 



(3.23) 



(3.24) 



we can write, for x > 0, 



Pr{X 3 (5i(t)) e cte}/d:r 



1 x 



q 37T y S 

2\fxl? 



Ki/a 



2 a; 3 / 2 \ t 1 /i 



3 3 / 2 / S 37T V s 



#1/3 



2 t 3 / 2 
(3.25) 



ds 



2Vxt* [°° 



2a; 3 / 2 
3 3 / 2 



s if 



1/3 



2t 3 ' 2 
W 2 



(3.26) 



In view of (see [2| formula 6.521]) 



sK v (ys) K v {zs) ds 



ir(yz) v (y 



2 sin7rz/ (y 2 — z 2 ) 



-, $l{y + z}>Q, \3l{v}\<l 



we get that 



3 1 / 2 t x-t 

Pr{X 3 (Si(t))edx} = dx — = ^ a;,t>0 (3.27) 

3 2tt r - r J 



which coincides with (3.18). 

The Cauchy densities pertaining to the composition X__x_(S __i_.(i)), t > 0, 
solve also a second-order p.d.e. as we show in the next theorem. 
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Figure 5. The profile of the function <|3.18 1. 
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Theorem 3.6. The Cauchy densities 

1 



f(x,t;m) = 



t COS 7T- 
2m 



7: ^ + ts[n ^L.)2 +t 2 C0S 2_^_ 



m e N 



satisfy the following second-order equation 
dt 2 



dx 2 



n . * d 2 f 

2 sin - — — — — , x £ 



2m dxdt 



t > 0. 



Proof. It is convenient to write (3.28) as a composed function 

1 u 

TT -7/2 I it )2 



where 
Since 



t cos ■ 



2m ' 



x + t sin 



2m 



cos 



df 

dt 

d 2 f 

dl = d£ 

dx dv 



^_dj_ 
2m du 

2m du 2 
and 



sm 



■K_dJ_ 

2m dv 

7T . 7T d 2 f 
2m 2m dudv 



^#7 

2m dv 2 



dx 2 dv 2 



and 



we have that 



dt 2 



dx 2 



8*1 
du 2 



2m du 2 dv 2 



dv 2 



„ • TT 7T d 2 f 

2 sm — cos — — — 
2m 2m dudv 



(3.28) 



(3.29) 



2m dv 2 



dv 2 
-2 sin 

-2 sin 



1 — cos 

7T d 



sm 



2m dv 

J^d_df_ 
2m dx dt 



2m 2m 

n df TT 3/ 

sin o - p~ + cos a~ 
2m dv 2m du 



„ . 7T 77 (9/ 

2 sin — cos 

2m 2m dudv 
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□ 



Remark 3.7. The characteristic function of (3.28) is 

e+OO 



J0X 



— t\j3\ cos ^——i/3t sin s^r 



f(x, t; m)dx = e 

and can be obtained by considering the bounded solution to the Fourier transform 

d 2 F „ 7T dF 



of (3.29) 



, , + 2iB sin - 
dt 2 P 2m dt 



0'F = 0. 

For the even-order Laplace equations we have the following result. 
Theorem 3.8. The solution to the higher-order Laplace-type equation 

x e K, t > 



d 2n u 



d 2n u 
dx 2n ' 



dt 2n 

subject to the initial conditions 
(u(x,0) = 5(x) 

is the classical Cauchy distribution given by 
u(x,t) = Pr{X 2n (SjL(t)) € dx}/dx = 



_ (-1)" fe! 7T(fc+l) 

fe+i COS 2 : 



t 



< k < 2n 



ir(x 2 + t 2 )' 

where X2 n (t), t > is a pseudo-process such that 

Ee i0X 2n (t) = R -tp\ 

Proof. The pseudo-process X2 n (t), t > related to the equation 

— = < iv t + 1 ^ 2 " u 

has fundamental solution whose Fourier transform reads 

e if3x u(x,t)dx = e' tp2n . 
If S_i_ (t), t > is a stable subordinator with Laplace transform 

Eexp (-AS_L(i)) = exp f-U^") , A > 0, £ > 
the characteristic function of X 2rl (S'_!_ (<)), t > becomes 



e ^P r {X 2n (5i(t))Gdx} 



Pr{Si(t)eds} 



(3.30) 



(3.31) 



(3.32) 



(3.33) 



= exp(-t\/3\e 1 ^) , r = 0,1,..., 2ti-1 (3.34) 

For r = 0, we have the characteristic function of the Cauchy symmetric law. For 
r =/= and n < r < 2n — 1 we have a function which is not absolutely integrable 
and, for 0<r<n— lis not a characteristic function (but can be regarded as a 
Cauchy r.v. at a complex time). The functions 

F r (p,t) = e~ t ^ ¥i 
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for all < r < 2n — 1 are solutions to 

d 2n F r 
dt 2n ~~ 



We now check that for < k < 2n — 1 the initial conditions (3.31) are verified by 
the Cauchy distribution. Indeed, 



dt k 



_ d k f 1 
t=o ~dj* 

^ r+oo 

_(~l) k k\ 
7r|x| fe+1 



+oo 



t=0 



e- i/3x (-l) k \(3\ k d/3 
ir(k + 1)' 



cos 



□ 



x e 



t > 



Remark 3.9. We notice that for n — 1 the problem above becomes 

d 2 u d 2 u 

W = "dx 2 ' 
subject to the initial conditions 

it(x, 0) = S(x) 



which is in accord with 



It ( x ' *) 
9 



tt|x| 



2" COS 7T 



t=0+ 



t 



dtn{x 2 +t 2 ) 



1 



t=o+ 7rx^ 



The connection between wave equations and the composition of two independent 
Cauchy processes C 1 (|C 2 (t)|), t > has been investigated in D'Ovidio and Ors- 
ingher [1 and more general results involving the Cauchy process have been pre- 
sented in Nane [8]. 



Remark 3.10. We finally notice that the equation 

d 6 u d 6 u 



dt e dx 6 



= 



can be decoupled as 



Q3 Q3 \ / Q3 Q3 



dt 3 dx 3 J V dt 3 dx 3 



u = 0. 



(3.35) 
(3.36) 



Form the Corollary 3.3 the solution to (3.36) can be therefore written as 

V3_4 



u(x, t) = 



1 

V3 

~2 2 7T 



V3 



te l s 



te e 
2 



+ \t 2 e-^ 



e if ( x 2 + t e -if + ^-if + |f2g-if j 



(x 2 + f e"*? + arte-*? + §t 2 e- l i) (x 2 + |e 4 i + xte lJ e + |f.Vf ) 

3+2, 



e l e (x 2 -I- |e*3 -I- xte l e + |i 2 e l 3 ) 



(x 2 + fe"^ + xie _i S + \t 2 e~ l l ) (x 2 + |e i: J + xie'f + ftVi) 
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^3^ (x 2 + t 2 ) cos ^ + xt 



2 2 7r (x 2 + te~^ + xte~ i ^+) (x 2 + te^ + xte i2 s) 
_V3 t {x 2 + t 2 )cos%+xt 

2 27r (x 2 + t 2 +xtcos^) 2 + 2x 2 t 2 cos^ ' 



Equation ( 3.36[ ) is satisfied by the Cauchy density and therefore by the probability 



law (3.371 which however is no longer a Cauchy distribution but is unimodal and 



asymmetric. 
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